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Microwave-radiation induced giant magnetoresistance oscillations recently discovered in high- 
mobility two-dimensional electron systems are analyzed theoretically. Multiphoton-assisted impurity 
scatterings are shown to be the primary origin of the oscillation. Based on a theory which considers 
the interaction of electrons with electromagnetic fields and the effect of the cyclotron resonance in 
Faraday geometry, we are able not only to reproduce the correct period, phase and the negative 
resistivity of the main oscillation, but also to predict the secondary peaks and additional max- 
ima and minima observed in the experiments. These peak-valley structures are identified to relate 
respectively to single-, double- and triple-photon processes. 
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The discovery of a new type of giant magnetorersis- 
tance oscillations in a high mobility two-dimensional 
(2D) electron gas (EG) subject to crossed microwave 
(MW) radiation field and a magnetic field ;^'^'?''*'^ es- 
pecially the observation of "zero-resistance" states de- 
veloped from the oscillation minima^*^*^ has revived 
tremendous interest in magneto-transport in 2D electron 
systems i2i2ii2iiiii2ii^ These radiation-induced oscillations 
of longitudinal resistivity Rxx are accurately periodical 
in inverse magnetic field 1 / B with period determined by 
the MW frequency u) rather than the electron density iVe. 
The observed Rxx oscillations exhibit a smooth magnetic- 
field variation with resistivity maxima at uj /uJc = i ~ <5- 
and minima at uj/lJc = j -f <5+ {oJc is the cyclotron fre- 
quency, j = 1,2,3...) having positive 5± around The 
resistivity minimum goes downward with increasing sam- 
ple mobility and/or increasing radiation intensity until a 
"zero-resistance" state shows up, while the Hall resistiv- 
ity keeps the classical form Rxy = B/Nf,e with no sign 
of quantum Hall plateau over the whole magnetic field 
range exhibiting Rxx oscillation. 

To explore the origin of the peculiar "zero- 
resistance" states, different mechanisms have been 
suggestediSiSiiLiLiSiii It is understood that the appear- 
ance of negative longitudinal resistivity or conductivity 
in a uniform model suffices to explain the observed van- 
ishing resistance)^ The possibility of absolute negative 
photoconductance in a 2DEG subject to a perpendic- 
ular magnetic field was first explored 30 years ago by 
Ryzhiiii^ii^ Recent worksSiiSiii indicated that the peri- 
odical structure of the density of states (DOS) of the 
2DEG in a magnetic field and the photon-excited elec- 
tron scatterings are the origin of the magnetoresistance 
oscillations. Durst et alA' presented a microscopic analy- 
sis for the conductivity assuming a 5-correlated disorder 
and a simple form of the 2D electron self-energy oscillat- 
ing with the magnetic field, obtaining the correct period, 
phase and the possible negative resistivity. Shi and Xie^^ 
reported a similar result using the Tien and Gorden cur- 
rent formulap^ for photon-assisted coherent tunneling. In 
these studies, however, the magnetic field is to provide an 



oscillatory DOS only and the high frequency (HP) field 
enters as if there is no magnetic field or with a magnetic 
field in Voigt configuration. The experimental setup re- 
quires to deal with the magnetic field B perpendicular to 
the HF electric field. In this Faraday configuration, the 
electron moving due to HF field, experiences a Lorentz 
force which gives rise to additional electron motion in the 
perpendicular direction. In the range of ~ uic, the elec- 
tron velocities in both directions are of the same order of 
magnitude and are resonantly enhanced. This cyclotron 
resonance (CR) of the HF current response will certainly 
change the way how the photons assist the electron scat- 
tering. 

In this Letter, we construct a microscopic model for 
the interaction of electrons with electromagnetic fields 
in Faraday geometry. The basic idea is that, under the 
influence of a spatially uniform HF electric field, the 
center-of-mass (CM) of the whole 2DEG in a magnetic 
field performs a cyclotron motion modulated by the HF 
field of frequency u. In an electron gas having impurity 
and/or phonon scatterings, there exist couplings between 
this CM motion and the relative motion of the 2D elec- 
trons. It is through these couplings that a spatially uni- 
form HF electric field affects the relative motion of elec- 
trons by opening additional channels for electron transi- 
tion between different states. Based on the theory for 
photon-assisted magnetotransport developed from this 
physical idea, we show that the main experimental results 
of the radiation-induced magnetoresistance oscillations 
can be well reproduced. We also obtain the secondary 
peaks and additional maxima and minima observed in 
the experimentsi^ 

For a general treatment, we consider electrons in a 
unit area of a quasi-2D system in the x-y plane with a 
confining potential V{z) in the z-direction. These elec- 
trons, besides interacting with each other, are scattered 
by random impurities/disorders and by phonons in the 
lattice. To include possible elliptically polarized MW il- 
lumination we assume that a uniform dc electric field Eq 
and ac field Et = Eg sin(a;t) -I- Ec cos{ujt) of frequency 
Lu are applied in the x-y plane, together with a mag- 
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netic field B = (0, 0, B) along the z direction. In terms 
of the 2D CM momentum and coordinate of the elec- 
tron system jiliiSii^ which are defined as P = J2j Pj\\ ^^d 
R- = ^e"^ J2j TjII with PjIi = {pjx,Pjy) and rj|| = {x.j,yj) 
being the momentum and coordinate of the jth electron 
in the 2D plane, and the relative electron momentum 



and coordinate p^|| = p,- 



-P/N^ and r' 
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R, the 



Hamiltonian of the system can be written as the sum of 
a CM part Hem and a relative- motion part Hd- (A(r) is 
the vector potential of the B field) , 

Hem = 77T^(P - iVeeA(R))2 - N^Eo + Et) • R, (1) 



^e.-E[i(p;i|-^^A(r;„) 



together with couplings of electrons to impurities and 
phonons, H^i and H^-p- Here m and rUz are respec- 
tively the electron effective mass parallel and perpendic- 
ular to the plane, and Vc stands for the electron-electron 
Coulomb interaction. Note that although in Hem and H^v 
the CM and the relative electron motion are completely 
separated, the CM coordinate R enters Hd and H^p^^^ 
Starting from the Heisenberg operator equations for the 
rate of change of the CM velocity \ ^ -i\V,H]+ dW/dt 
with V = — ?[R, H], and that of the relative electron en- 
ergy 7?er = H], we proceed with the determina- 
tion of their statistical averages. 

The CM coordinate R and velocity V in these equa- 
tions can be treated classically, i.e. as the time-dependent 
expectation values of the CM coordinate and velocity^^ 
R(t) and V(t), such that R(t) - R(t') = V(s)ds. We 
are concerned with the steady transport under an irradia- 
tion of single frequency and focus on the photon-induced 
dc resistivity and the energy absorption of the HF field. 
These quantities are directly related to the time-averaged 
and/or base-frequency oscillating components of the CM 
velocity. At the same time, in an ordinary semiconductor 
the effect of higher harmonic current is safely negligible 
for the HF field intensity in the experiments. Hence, it 
suffices to assume that the CM velocity, i.e. the electron 
drift velocity, consists of a dc part vq and a stationary 
time-dependent part that 



V(i) Vq + vi cos{ijjt) + V2 sin(cji). 



(3) 



This time-dependent CM velocity enters all the operator 
equations having couplings to impurities and/or phonons 
in the form of the following exponential factor, which can 
be expanded in terms of Bessel functions J,i(x): 

n— — oo 
oo 



Here ^ = ■ Vi)^ -|- (q|| • \2Y/^ and ianip = (q • 

V2)/(q • vi). On the other hand, for 2D systems hav- 
ing electron sheet density of order of 10^^ m^^, the 
intra-band and inter-band Coulomb interactions are suffi- 
ciently strong that it is adequate to describe the relative- 
electron transport state using a single electron temper- 
ature Te. Except this, the electron-electron interaction 
is treated only in a mean-field level under random phase 
approximation (RPA)iiSiiS For the determination of un- 
known parameters Vq, Vi, V2, and Te, it suffices to know 
the damping force up to the base frequency oscillating 
term F(i) — Fo-l-Fs sin{ijot) + ¥ c cos(wi), and the energy- 
related quantities up to the time-average term. We finally 
obtain the following force and energy balance equations: 
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N^eEo -vo + Sp-W = 0. 



(4) 
(5) 

(6) 
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Here 



Fo = E 1^^(^11)1' E qi|J'(On2(q||,c.o-nc.) 

qy n— — oo 

oc 

+ El^^(q)l' E qi|J'(OA2(q,^o + f^q-nc.) (8) 



n— — oo 



is the time-averaged damping force, Sp is the time- 
averaged rate of the electron energy-gain from the HF 
field, ^Nce{Es • V2 -I- Ec ■ Vi), which can be written in 
a form obtained from the right hand side of Eq. ((SJ by 
replacing the qy factor with nuj, and W is the time- 
averaged rate of the electron energy-loss due to coupling 
with phonons, whose expression can be obtained from the 
second term on the right hand side of Eq. ((SJ by replac- 
ing the q|| factor with Qq, the energy of a wavevector- 
q phonon. The oscillating frictional force amplitudes 
F^ = F22-F11 and F^ = F2i-fFi2 are given by (/z = 1,2) 

00 

F,^^-J2m\ui<in)\' E [4(0]'ni(q||,^o-nc^) 

qy n— — 00 

00 

-Eqil^Ml^^(q)P E K(0]'Ai(q,u;o + r!q-na;), 



n— — 00 



F2^ = Eq|l7l^(q||)l' ^ 2nj2(OH2(q||,^o-n^) 

qii 



n— — 00 



n— — 00 



+ ^q|,^|M(q)p J2 2nj2(^)A2(q,c.o + r!q-7ic.). 

q n— — OQ 

In these expressions, ry^^ = qy ■ v^/uj£^; cjq = qy • 
Vq; C/(qy) and M{q) stand for effective impurity and 
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phonon scattering potentials, n2(q||,f2) and A2(q, f2) = 
2n2(q||,rJ)[n(riq/r)-n(r2/re)] (withn(x) = l/(e^-l)) 
are the imaginary parts of the electron density correla- 
tion function and electron-phonon correlation function 
in the presence of the magnetic field. ni(q||,i7) and 
Ai(q, il) are the real parts of these two correlation func- 
tions. The effect of interparticle Coulomb interactions 
are included in them to the degree of level broadening 
and RPA screening. 

The HF field enters through the argument ^ of the 
Bessel functions in Fg, F^i,, W and Sp. Compared with 
that without the HF field {n = term on\y)^ we see 
that in an electron gas having impurity and/or phonon 
scattering (otherwise homogeneous), a HF field of fre- 
quency io opens additional channels for electron transi- 
tion: an electron in a state can absorb or emit one or 
several photons and scattered to a different state with 
the help of impurities and/or phonons. The sum over 
n| > 1 represents contributions of single and multiple 
photon processes of frequency-w photons. These photon- 
assisted scatterings help to transfer energy from the HF 
field to the electron system (S'p) and give rise to addi- 
tional damping force on the moving electrons. Note that 
Vi and V2 always exhibit CR in the range w ~ Wc, as can 
be seen from Eqs.© and © rewritten in the form 



vi = (i-.,v-r^{— 

L muj 



1 



N^muj 



V2 = {u;'Ju;' - I)-' { 
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(10) 



Therefore, ^ may be significantly different from the ar- 
gument of the corresponding Bessel functions in the case 
without a magnetic field or with a magnetic field in Voigt 
configuration.^'^ 

Eqs. Q-lQ can be used to describe the transport and 
optical properties of magnetically-biased quasi-2D semi- 
conductors subject to a dc field and a HF field. Taking 
vo = {vox,0,0) in the x direction, Eq. Q yields trans- 
verse resistivity Rxy = Eoy/NeCVox = B/N^e, and lon- 
gitudinal resistivity Rxx = Eox/Ngevox = -Fo/N^e'^vox- 
The linear magnetoresistivity is then 



rtxx - 2^qx ^2^2 ''nK^) 



|M(q)| 



n— — oo 
2 oo 



(11) 



£1— H-noj 



The parameters vi, V2 and To in (|ll|l should be deter- 
mined by solving equations Q, © and (0) with vanish- 
ing vq. We see that although the transverse resistivity 
Rxy remains the classical form, the longitudinal resistiv- 
ity Rxx can be strongly affected by the irradiation. 



We calculate the unscreened H2(q||, ft) function of the 
2D system in a magnetic field by means of Landau 
representation^ii 

n2(q||,0) = J—J2Cn.n'{lBq\\/m2{n,n,n), (12) 



Il2{n,n',n) 



de [/(e) - /(£ + n)] 



X ImG„(e-|-r2)ImG'„'(£), (13) 



where I-q = yJl/\eB\ is the magnetic length, C„_„+i(F) = 
n\[{n+iy]-'^Y^ e^^ [L^„{Y)Y withij,(y) the associate La- 
guerre polynomial, /(e) = {exp[(e — m)/^o] + 1}^^ the 
Fermi distribution function, and ImG'„(e) is the imagi- 
nary part of the Green's function, or the DOS, of the 
Landau level n. The real part functions Hi(q||,17) and 
Ai(q|| , J7) can be obtained from their imaginary parts via 
the Kramers-Kronig relation. 

In principle, to obtain the Green's function G'„(e), 
a self-consistent calculation has to be carried out with 
all the scattering mechanisms includediSi In this Letter 
we do not attempt a self-consistent calculation of G„(e) 
but choose a Gaussian-type function for the purpose of 
demonstrating the Rxx oscillations (£„ is the energy of 
the n-th Landau level) 

ImG„(c) - -(7r/2)^r-i exp[-(e - e„)V(2r2)] (W) 

with a broadening parameter F (2eci;ca/7i''7i/io)^^^, 
where is the linear mobility at temperature T in the 
absence of the magnetic field and a > 1 is a semiem- 
pirical parameter to take account the difference of the 
transport scattering time determining the mobility /io, 
from the single particle lifetimeASiiS 

The moderate microwave intensity for the Rxx oscilla- 
tion in these high-mobility samples yield only slight elec- 
tron heating, which is unimportant as far as the main 
phenomenon is concerned and is neglected for simplic- 
ity. We consider scatterings from remote impurities as 
well as from acoustic phonons. After solving vi and V2 
from Eqs. jnj and Ij 10(1 the magnetoresisivity Rxx can be 
obtained directly from Eq. Hll|l . At lattice temperature 
T = 1 K, the contribution from photon-assisted phonon 
scattering is minor. The role of acoustic phonons, how- 
ever, becomes essential at elevated lattice temperatures. 
Calculations were carried out for linearly polarized MW 
fields with multiphoton processes included. 

Fig. 1 shows the calculated longitudinal resistivity Rxx 
versus uj/uJc = 7c subject to a linearly polarized MW 
radiation of frequency uj/2Tr — 100 GHz at four values 
of amplitude: i?s = 20, 45, 65 and 80 V/cm. Shubnikov- 
de Haas (SdH) oscillations show up strongly at high ojc 
side, and gradually decay away as 1/ujc increases. All 
resistivity curves exhibit pronounced oscillation having 
main oscillation period 7c = 1 (they are crossing at inte- 
ger points 7c = 2,3,4,5). The resistivity maxima locate 
around 7c = j — and minima around 7c = j + with 
5± ~ 0.23 - 0.25 for j = 3,4,5, S± - 0.17 - 0.21 for 
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FIG. 1: The longitudinal magnetoresistivity R^x of a 
GaAs-based 2DEG subject to a lineraly polarized HF field 
EaSiniuii). The parameters are: temperature T = IK, elec- 
tron density = 3.0 x 10^^ cm"" 
mobility = 2.4 x 10^ cm^ V" 
coefficient a = f 2. 



zero-magnetic-field linear 
s~^, and the broadening 



J = 2, and even smaller 6± for j — 1. The magnitude 
of the oscillation increases with increasing HF field in- 
tensity for 7c > 1.5. Resistivity gets into negative value 
for Eg = 80 V/cm around the minima at j — 1,2 and 
3, for Es = 65 V/cm at j = 1 and 2, and for Es = 20 
and 45 V/cm at j = 1. All these features, which are in 
fairly good agreement with experimental findings^*^*^ 
are relevant mainly to single-photon (|n| — 1) processes. 
Anomaly appears in the vicinity of 7c = 1, where the 
CR greatly enhances the effective amplitude of the HF 
field in photon-assisted scatterings and multiphoton pro- 



cesses show up. The amplitudes of the j = 1 maximum 
and minimum no longer monotonically change with field 
intensity. Furthermore, there appears a shoulder around 
7c — 1.5 on the curves of Eg = 45 and 65 V/cm, and 
it develops into a secondary peak in the Eg — 80 V/cm 
case. This has already been seen in the experiment (Fig. 2 
in Ref.0). The valley between 7c = 1.4 and 1.8 peaks 
can descend down to negative as Es increasing further 
(Fig. lb). The appearance of the secondary peak is due 
to two-photon (|n| = 2) processes. 

Radiation-induced resistivity behavior at 7c < 1 is 
shown more clearly in the cj/27r = 60 GHz case. As 
seen in Fig. Ic, a shoulder around 7c = 0.4-0.6 with a 
minmum at 7c — 0.6 can be indentified from the SdH os- 
cillation background for all three curves, which is related 
to two-photon process. With increasing MW strength 
there appears a clear peak around 7c = 0.68 and a val- 
ley around 7c = 0.76. This peak-valley is mainly due to 
three-photon (|n| = 3) process. In the case of 40 GHz, 
similar peak and valley also show up (Fig. Id). 

Qualitatively, the main oscillation features come from 
the symmetrical property of the DOS function in a mag- 
netic field. Since G„(£— jWc) = Gn-jis) for any integer j, 
the impurity contribution to Rxx from the rt-photon pro- 
cess, which is related to the weighted summation of the 
derivative H2 function over all the Landau levels at fre- 
quency nuj [Eq. (|ll|l ] , has an intrinsic periodicity charac- 
terized by nu — jujc- The main oscillation of Rxx shown 
in Fig. la relates to single-photon process and character- 
ized by a; = jojc- We have also performed calculation us- 
ing a Lorentz-type DOS function and find that, although 
the oscillating amplitude and the exact peak and val- 
ley positions are somewhat different, the basic periodic 
feature of the radiation-induced magnetoresistivity oscil- 
lation remains. 
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